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We study charge ordered d-wave resonating valence bond states (dRVB) in the doped cuprates, 
and estimate the energies of these states in a generahzed t — J model by using a renormalized mean 
field theory. The long range Coulomb potential tends to modulate the charge density in favor of 
the charge ordered RVB state. The possible relevance to the recently observed 4x4 checkerboard 
patterns in tunnelling conductance in high Tc cuprates is discussed. 
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A number of recent STM experiments have shown spa- 
tial modulations in tunnelling conductance in high Tc 
cuprates More recent low temperature STM 

experiments have reported bias independent modulations 
of period approximately 4 — 4.5a (a: lattice constant) in 
the tunnelling conductance over a wide range of energy on 
underdoped 5^2212 ,6] and NaCCOC |7|]. Several theo- 
retical proposals have been put forward to interpret the 
observed checkerboard charge ordering |^ 0, ITol IllL IT^ . 
Chen et al ^ have proposed that the modulations are 
due to Cooper pair density wave. Fu et al have exam- 
ined the possibility of a soliton crystal in a generalized 
Hubbard model including a nearest neighbor Coulomb 
repulsion and an antiferromagnetic spin exchange cou- 
pling. Anderson p^] has proposed an explicit wavefunc- 
tion describing a Wigner solid of holes embedded in a 
sea of (i-RVB state, and pointed out that the long-range 
Coulomb interaction furnishes the energy gain and the 
stiffness of the hole wavefunction opposes the deforma- 
tion. The detailed calculations, however, have not been 
carried out in Ref. [T^ . 

In the present paper, we study the charge ordered 
dRVB in the doped cuprates. We use a Gutz wilier pro- 
jected wavefunction with both BCS pairing and charge 
ordering to describe the charge ordered state in cuprates. 
Our approach is similar to the idea outlined by Ander- 
son who formulated the charge ordering in a dRVB 
by a site-dependent fugacity, which was introduced by 
Laughlin in the context of Gutzwiller projected state in 
study of the Gossamer superconductivity 0, . Here 
we shall use the renormalized mean field theory (RMFT) 
developed early 0, to formulate charge ordering by 
site-dependent renormalization factors and estimate the 
energies of these states in the t — J model. We show that 
the long range Coulomb potential tends to modulate the 
charge density in favor of the charge ordered RVB state 
and that the favorable patterns for the charge ordering 
depend on the doping concentration. Our calculations 
suggest that the observed checkerboard patterns may well 
be induced by the long range Coulomb repulsion. 

We consider a generalized t — J model with an addi- 



tional long-range electron Coulomb potential. 



H = 

Ht = 

Hs = 

Hr = 
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where Qct is an annihilation operator of a spin a elec- 
tron at site i. The sums in the kinetic and spin-exchange 
Hamiltonian Ht and Hs run over all the nearest neighbor 
pairs, rii = X^^n^^, and n^^ = c\^Cia. The sum in He 
runs over all the sites of i and j. e is the dielectric con- 
stant, and Vij is the spatial distance between the two sites 
i and j. An positive charge background to balance the 
charge neutrality is implied. There is a local constraint 
on every site, cJ^Qcr < 1. In this Hamiltonian, He fa- 
vors a charge ordering, while Ht prefers a uniform charge 
distribution. 

We use the RMFT to estimate the energies of the 
charge ordered RVB states. The RMFT was developed 
for the t — J model to study a charge homogeneous RVB 
state 0, • ^^^^ theory, one considers a Gutzwiller 
projected BCS state for the possible superconducting 
(SC) ground state. Here we shall extend it to the charge 
inhomogeneous case. We consider a Gutzwiller projected 
state, 



|^)=Pg|^o), 



(1) 



where | ^o) is a charge ordered BCS state, and Pq = 
H^(l — rii^riii) is the Gutzwiller projection operator. We 
use Gutzwiller 's approximation to relate the expectation 
values of the kinetic or spin exchange energies in the pro- 
jected state I ^) (denoted by <>) to the corresponding 
expectation values in the unprojected state | ^o) (de- 
noted by <>o) by two different renormalization factors 
Qt and Qs'. 



(2) 
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The renormalization factors are determined by the ratio 
of the probabihties of the physical process in the pro- 
jected and in the unprojected states [T^, [la IlS Il3 • Sim- 
ilar to the method used for the homogeneous case, 
we find 

ij ^ ^ / (l-nO(l-T^ 

.. 4 

" (2-n,)(2-n,)* 

They depend on the electron densities at the sites i and 
j. In the homogeneous case, Ui = g' s are independent 
of the sites, we recover the results in Ref. — 
2x/(l + x) and gf^ = 4/(1 + x)^, with x = I — n the 
hole density. The variational calculation of the projected 
state I ^) in is then mapped onto the unprojected state 
I ^o) in a renormalized Hamiltonian i^e//? given by 

H[ = -tY^ gl^ cl^Cja ^ h.c. 

H's = J^gl^S.-S,. (4) 

(id) 

Note that the intersite Coulomb interaction is not renor- 
malized in the theory {i.e. the renormalization factor is 

!)■ _ 

Similar to the procedure in Ref. we introduce 

two mean fields: a particle- hole amplitude field ^ij = 

5^cr(4cr^jcr)o5 a particle-particle pairing field Aij = 
{ci^Cji — CiiCj^)o. The renormalized Hamiltonian can 
then be solved by a self-consistent mean field theory. The 
energy of Heff in the unprojected state, hence the en- 
ergy of the generalized t — J model in the projected state 
can be written in terms of the self-consistent mean fields. 

In the uniformly charged dRVB state, = ^, and Aij = 

±A. The energy per site is found to be E^^'' = —Atg^.^'^^,— 

(3J/4)^s^^(^^ + A^), where we have dropped out the long 
range Coulomb energy of a uniform electron density for 
it cancels to the energy due to the oppositely charged 
background. 

In the inhomogeneous case, the self-consistent equa- 
tions, or the Bogoliubov de-Genes equations, are more 
complicated. In what follows we shall make an approx- 
imation to replace the mean fields and Aij by their 
average mean values obtained in the uniform dRVB state, 
and consider the effect of charge ordering on the kinetic 
and spin-exchange energies due to the renormalization 
factors gl^ and gl^ ^ and on the Coulomb potential. This 
is a rather dramatic approximation similar to what pro- 
posed by Anderson [13 , but it should capture a substan- 
tial part of the effect of the charge ordering. The accu- 
racy of this approximation will be examined in a limiting 
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FIG. 1: The mean fields A and ^ versus the hole concentration 
X in the uniform dRVB state. t/J = 3. 

case, which turns out to be quite good. Within this ap- 
proximation, the energy per site of the charge ordered 
dRVB state relative to the uniform dRVB state is, 

AE = AEt + A^;^ + AEc 

AEt = {gt-g't){Ht)o/Ns 

AEs = {gs-g^s{Hs)o/Ns 

AE, = eV(2e7V,)^(n,n, -n')/r,,. (6) 

In the above equations, gt^s = 9t%/^^s, {Ht,s)o is 

the average kinetic (spin exchange) energy in the uniform 
dRVB state. In practice, we first solve the RMFT for the 
uniform dRVB state, from which we obtain ^, A, and 
{Ht^s)o- In Figffl we plot the two mean fields as func- 
tions of hole doping x for J/t = 1/3. We then calculate 
gt^s and AE, for various types of charge ordering patterns 
to estimate the energy of the charge ordered RVB state, 
and to determine the optimal charge distribution. The 
calculation of the long range Coulomb energy is similar 
to the calculation of the Modelung constant, and it con- 
verges rapidly with the appropriate choice of the method 
in summation. 

Motivated by the approximate 4x4 charge ordered 
states observed in STM experiments, we consider several 
types of parent patterns shown below with periodicity 
of 4a along both directions in the square lattice. Each 
symbol represents a lattice site, and the sites marked with 
the same symbol have the same electron density. 
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We denote = n + ^1, Uo = n -\- 62^ and = n + ^3. 
Since the overall average electron density of the system is 
n, only two out of the three S^s are independent. By using 
Eq. (4), we have gl, = + ^^M' C = + JdZ' 
S't'i = bn + lan + hn+hVs, and^t^ = lg** + lg°° + 
+ I^M + i^M' ^^^^ the superscript in g indicates 
the type of the parent pattern, and the superscript in 
g refers to the two sites with the marked symbols. The 
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FIG. 2: Low energy charge ordering 
pattern Al of 4a x 4a symmetry at hole 
density x around 1/16. Shown are 9 x 
9 patches. Each square represents a 
lattice site and the light (dark) square 
represnts low (high) electron density. 



Coulomb energy can be shown to be quadratic ind^s, and 
they are given by, in unit of e^/ea, 



TABLE L The energy and charge distribution of pattern Al at 
hole density x = 0.0625 and x = 0.05. An ideal hole crystal state 
with riB = 1 at a: = 0.0625 is also listed for comparison. 



= -4.039(5? - 4.039(5i - 4.M7SiS2 
- -1A9751 - 7.95951 - 3.468(51^2 
A£;f^ = -0.567(5? - 3.772(5^ - 2.5II6162 
AE^ = -1.248(5? - 4.138^1 - L428^i^2. (7) 

Using these expressions, we have optimized the energy 
by varying parameters Si and 62 , and obtained charge or- 
dered states with lower energies. These states are deriva- 
tives of the parent patterns under the consideration, but 
may have a higher symmetry than the parent state for 
those sites marked with different symbols may have the 
same electron density. Below we shall discuss our results 
in three different regions of the hole concentration. In 
all of our calculations, we use J/t = 1/3, t = 0.3eF, 
a = 3.8A. 

At the hole density around 1/16, the lower energy 
charge ordering pattern is Al as shown in FigO There 
are only two types of the distinct sites in terms of the 
electron density in this pattern. The numerical values of 
the energy gain and the charge distributions are given 
in Table ID for x = 1/16 and x = 0.05. Ah other pat- 
terns at these dopings have energies either higher than 
or too close to the energy of the uniform dRVB state 
{AE > -O.OleV), and are not listed here. At x = 1/16, 
and e = 1, the lowest energy state has a charge distri- 
bution slightly deviated from a commensurate state with 
the light site completely empty (n = 0) and the dark 
site fully occupied (n = 1). At a larger e, the energy 
gain decreases rapidly and the energy of the pattern Al 
is just slightly lower than that of the homogenous case 
at e = 1.5. There is no stable charge ordering pattern 
at e much larger than 1.5. x = 1/16 is an ideal hole 
density for the pattern Al, which was also discussed in 
im and suggested in the magnetic and optical measure- 
ments 0,123 • At X = 0.05, the pattern Al is stable only 
for smaller e, but is no longer stable for e = 1.5. Note 
that the pattern Al at x = 0.05 is an insulator for there 
is no any connected path for holes to move through the 
lattice. 

At the hole density around 1/8, there are several charge 
ordering patterns as shown in Fig|31 Among them the 
favorable pattern is 51. Patterns B2 and have three 
types of distinct sites in terms of the electron density, 
while Bl and B3 have two types of distinct sites. The 
energy gain and the charge distribution are given in Table 
^ for X = 1/8 and x = 0.1. Here we only list those 
patterns with relatively lower energies. As we can see 
from table in at x = 1/8 the energies of patterns Bl and 
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FIG. 3: Lower energy charge ordering patterns for hole density 
around x = 1/8. Bl is of a symmetry of s/Sa x \/8a, and B3 is a 
stripe. 



B2 are slightly lower than that of the homogeneous case 
ate = 2. Atx = 0.1, the energy gain due to the charge 
ordering at e = 2 is already very tiny. 

It is interesting to note that around the low hole den- 
sity X = 1/8, both the checkerboard pattern B2 and the 
stripe pattern|2l| B3 are SC states for holes in these pat- 
terns can move through the lattice. 

At high hole concentrations, several new charge order- 
ing patterns with lower energies appear, which are shown 
in Fig^ In Table IIIIl we list the energies and charge 
densities of the lower energy patterns at x = 0.15. For 
e = 1, the five patterns (51, (71, (72, (73, and (74) have 
very close energies. In the pattern C's, the electron den- 
sities at the dark and grey sites are quite close. The 
empty sites in patterns Bl and C's form a \/8a x ^/Sa 
Wigner hole crystal. Among the series C, pattern CI is 
a conducting phase. We do not find any lower energy 
charge ordering pattern ate>2.5. Atx = 0.2, the most 
favroable patterns are (71, (74 and the stripe pattern B3. 

In the energy estimation for the charge ordered RVB 
state, we have focused on the effect of the charge den- 
sity dependent renormalization factors, but neglected the 
site-dependence of the mean fields ^ and A. This ap- 
proximation turns out to be quite good in a limiting case 
where all the holes are completely localized at a single 
site. Consider pattern Al at x = 1/16 and pattern Bl 
at X = 1/8 with the electron density n either zero or 
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FIG. 4: Several lower energy charge ordering patterns at x = 0.15, 
not included in Figs 3. Patterns CI and C2 are related by the 
interchange of the dark and gray sites, so for patterns C3 and C4. 
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TABLE II: The energy and charge distribution of lower energy charge ordering patterns at hole density x = 0.125 and x = 0.1. 
The energies of some patterns with charge distribution of na = 1 are also listed for comparison. 
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TABLE III: Energies and charge distributions of lower energy charge ordering patterns at x = 0.15 
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1. In this limit, the kinetic energy vanishes. The spin 
exchange energy of the state can be estimated by a di- 
rect counting of the missing bonds due to the vacancies 
in an otherwise half filled background, which is given by 
Es = 2(1 — 2x)a per site, with a = — 0.344J the spin ex- 
change energy per bond at the half filling. For J = O.leV 
we have Es = -0.060eV at x = 1/16 and Es = -0.052eV 
at X = 1/8, which are very close to the results obtained 
in the present MET: Es = OmSeV at x = 1/16 and 
Es = -0.054eV at x = 1/8. 

In summary we have studied the charge ordered RVB 
states in the doped cuprates within a generalized t — J 
model by using a renormalized mean field theory. While 
the kinetic energy favors a uniform charge distribution, 
the long range Coulomb repulsion tends to spatially mod- 
ulate the charge density in favor of charge ordered RVB 
states. Since both the Coulomb potential and the leading 
order in kinetic energy are quadratic in the density vari- 



ation, we expect and indeed have found that the charge 
density variation from the uniform state is always large 
in the charge ordered state. The stability of the charge 
ordered RVB state strongly depends on the dielectric con- 
stant e. In cuprates, e ^ 2.5 — 5. Our calculation suggests 
that the observed charge ordered state in STM exper- 
iments in cuprates may well be interpreted due to the 
long range Coulomb interaction. Among the favorable 
charge ordered superconducting states, pattern Bl has a 
symmetry of \/8 x >/8, patterns B2 and CI both have 
checkerboard structure, and pattern B3 is a stripe. Be- 
cause of the intersite Coulomb repulsion, we do not find 
the bound hole paris in the charged ordered states. 
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